Electron Exchange Coupling for Single Donor Solid-State Spin Qubits 
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Inter-valley interference between degenerate conduction band minima has been shown to lead 
to oscillations in the exchange energy between neighbouring phosphorus donor electron states in 
silicon 1 '-. These same effects lead to an extreme sensitivity of the exchange energy on the relative 
orientation of the donor atoms, an issue of crucial importance in the construction silicon-based spin 
quantum computers. In this article we calculate the donor electron exchange coupling as a function of 
Oh! donor position incorporating the full Bloch structure of the Kohn-Luttinger electron wavefunctions. 

It is found that due to the rapidly oscillating nature of the terms they produce, the periodic part 
of the Bloch functions can be safely ignored in the Heitler-London integrals as was done by Koiller 
QQ ■ e * a! 1 ' 2 , significantly reducing the complexity of calculations. We address issues of fabrication and 

calculate the expected exchange coupling between neighbouring donors that have been implanted 
into the silicon substrate using an 15keV ion beam in the so-called 'top down' fabrication scheme 
for a Kane solid-state quantum computer. In addition we calculate the exchange coupling as a 
^ ■ function of the voltage bias on control gates used to manipulate the electron wavefunctions and 

t ' implement quantum logic operations in the Kane proposal, and find that these gate biases can be 

used to both increase and decrease the magnitude of the exchange coupling between neighbouring 
donor electrons. The zero-bias results reconfirm those previously obtained by Koiller et a&. 
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Solid state systems have emerged as a promising candidate for the construction of a large scale quantum computer 
(QC). In particular spin based architectures take advantage of the relatively long spin dephasing times of donor 

£h ! electrons or nuclei in silicon. Single qubit operations are performed by tuning the spin Zeeman energy to be resonant 
with an oscillating field which drives the transition while neighbouring qubits are coupled via the electron exchange 
interaction, whether it be directly in the case of electron-spin proposals^, or indirectly in the case of nuclear spin 

^ • quantum computers. 

In this work we concentrate on the Kane^ concept of single phosphorus donor nuclear spin qubits embedded in a 
silicon substrate, which is a leading candidate in the search for a scalable QC architecture. The Kane architecture 
(shown in Fig. [2) calls for the placement of phosphorus donors at subsitutional (face-centred cubic (fee)) sites in the 
host silicon matrix, with inter-donor spacings of order 200A. Quantum logic operations on the nuclear-spin qubits 
arc implemented through coherent control of the donor electron wavefunctions which are coupled to the donor nuclei 
through the contact hyperfine interaction. This control of the electron wavefunctions is achieved through application 
of voltage biases to control gates placed on the substrate surface above (A-gate), and between (J-gate) the phosphorus 
donors, which create electrostatic potentials within the device thus altering the form of the electron wavefunctions. 

In a recent paper, Koiller, Hu and Das Sarma (KHD)i presented theoretical evidence for oscillations in the electron- 
mediated exchange coupling as a function of inter-donor separation, and a strong dependence of this coupling on the 
relative orientation of the two donors with respect the structure of the silicon substrate. They calculate the exchange 
coupling using an approximate version of the Heitler-London formalism and by essentially ignoring the periodic part 
of the Bloch wavefunctions in the expression for the donor electron wavefunctions. In this article we eliminate both 
these approximations, by calculating the exchange coupling in the full Heitler-London formalism and by including the 
full Bloch structure of the donor electron wavefunction. We show that while the first approximation breaks down for 
small donor separations, as discussed by KHD themselves^, the second approximation, that of ignoring the periodic 
part of the Bloch functions, is in excellent agreement with the full calculations regardless of the donor orientation. 

This paper is organised as follows. In section [H] we review the calculation of the valence electron wavefunctions for 
phosphorus donors in silicon, in the Kohn-Luttinger effective mass formalism. In section ITTT1 we discuss the Heitler- 
London formalism used to calculate the exchange coupling between neighbouring donor electrons, and we show that 
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FIG. 1: The Kane architecture based on buried phosphorus dopants in a silicon substrate 



the terms that arise due to the periodic part of the Bloch functions oscillate sufficiently rapidly so as to average to 
zero over the range of the integrals. Section IIVI contains a discussion of some of the fabricational issues that arise 
from the strong position dependence of the exchange energy on the production of a phosphorus donor based solid 
state quantum computer. 

Section [3 is devoted to a calculation of the dependence of the exchange coupling between neighbouring donor 
electrons on voltage biases applied to the control J-gate, used to tune the inter-qubit coupling. By extending the 
conventional effective mass formalism we construct a basis of non-isotropic hydrogen like envelope functions in which 
we expand the donor-electron wavefunction in the presence of the electric potentials. The potential created inside 
the device due to the J-gate bias is calculated using a commercial software package specifically designed for the 
modelling of semiconductor devices. Donor wavefunctions are obtained by direct diagonalisation, and the Hcitlcr- 
London formalism is used to determine the exchange coupling for various gate biases and donor separations oriented 
in the silicon [100] direction. 



II. THE DONOR ELECTRON WAVE FUNCTION 



Although the qubits of the Kane quantum computer are encoded by the nuclear spins, it is the donor electron 
which mediates both single and coupled gate operations. Therefore, the crucial element in the quantum description 
of the device is the donor electron wave function. Initial attempts to describe the operation of the device, particularly 
in response to external time dependent gate potentials which complicate the situation considerably, have focused on 
effective hydrogenic-type approximations for the donor electron wave function^, or the use of simplified potentials 7 . 
These calculations provide some useful estimates, however more detailed calculations are required, using realistic 
models of both the donor electron wavefunction and the potentials induced inside the device by the application of 
control gate biases. 

In going beyond the effective hydrogenic treatment of the ground-state electron wave function for phosphorus donors 
in silicon it is recognised that the underlying crystal structure of the silicon must have an effect. The wavefunction is 
thus expanded in the basis of the Bloch functions for siliconj&S 

*(r) = J F(k)0 k (r)dk. (1) 

The coefficients F(k) are obtained by substituting the above form into the Schrodinger equation, with the Hamiltonian 
H = Ho — U(r), where Hg is the Hamiltonian for the pure silicon crystal and U(r) is the donor-potential for phosphorus. 
The Bloch functions can be written in the form 4>k( r ) — e «k( r )> where itk(r) is a function that shares the periodicity 
of the lattice, and can be expanded in a basis of plane waves with wave-vectors equal to the reciprocal lattice vectors 
for silicon G, Uk(r) = ^ G Ak(G)e tGr . The result is that the Schrodinger equation can be written as 

EF(k) = E^F(k) +J2f K>,G>Ak,GU(r)e l ^- ] <'^ r e l( - G - G '') r F(k , )drdk' 

G,G' •* 

= E°F(k) +Y, I ^, iG ,A KG U{k + G — k' — G')F(k')dk', (2) 

G.G' J 
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where U(k) = J t/(r)e _lr k dr is the Fourier transform of the impurity potential. The E^. are the eigenenergies of the 
Bloch functions </>k( r ), for the pure silicon lattice. We make the approximation that due to the increased energies of the 
higher bands, only conduction band states contribute to the impurity wavefunction. Further, in silicon there are six 
degenerate conduction band minima, located along the (100) directions in k-space, 85% of the way between the centre 
(r point), and the zone boundary (X point). Because of the reduced energies in these regions the envelope functions 
can be expressed as a sum of functions localised around each of the conduction band minima F(k) = -^mOO- 
In the effective-mass treatment 8,9,10 the Bloch energies are expanded to second order around the conduction band 

minima E^ s» ^- I + J- J , where kj^ is the component of k — k M perpendicular to k M and k\\ is the parallel 

component. The m±,mu are effective masses and the inequality of these two values reflects the anisotropy of the 
conduction band minima. An additional approximation is made whereby only the terms with G = G' in the potential 
term of the Schrodinger equation are included. The assumption is that U(k — k' + G — G') << U(k — k') for G ^ G'. 
This approximation is well satisfied for a Coulombic potential U(k) = l/(K7r 2 |k| 2 ), with k — 11.9 the dielectric 
constant for silicon, and with the reciprocal lattice vectors of magnitude |G| = n2ir/d, where n is an integer and 
d = 5.43A is the lattice spacing for silicon. Another way of viewing the approximation is that the terms e^ G_G ' r 
in the first line of Eq. [5] oscillate on a scale sufficiently short compared to variations in U(r), that they average the 
integrand to zero. We will show in the next section that the same approximation allows us to ignore the periodic part 
of the Bloch functions in the donor electron wavefunctions when performing Heitler-London integrations. 

One further approximation is necessary, g ss ^4k,G, which coupled with the relation ^ G |^4k,G| 2 = 1 gives the 
effective mass Schrodinger equation 



% + ^) E w + J^ k - k ') E F ^ dk ' = ^E^(k). 



(3) 



In the standard effective mass treatment the so-called valley-orbit terms, which couple the envelope functions at 
different conduction band minima are ignored, and we are left with six independent equations, one for each minimum. 
With a Coulombic impurity potential the solutions are non-isotropic hydrogenic wavefunctions of the form 

F± M = exp[-y/(x 2 +y 2 )/a± + z*/a\\] ^ 

where F^(r) = /F M (k — k Al )e lk,r dk, and we have used the example of the envelope function localised around the 
2-minima of the conduction band. The Kohn-Luttingei4i& form for the electron wavefunction of a donor situated at 
any position R is then given by 

*(r - R) = J2 ^(r - R)e k -( r - R ^(r), , (5) 

where the periodic part of the Bloch function is indepedant of the position of the substitutional impurity. The values 
a±_ = 25.09A, a|| = 14.43A, are the effective Bohr radii, and are determined variationallyi^. 

It is well known that the valley orbit coupling does contribute to the the energy of the stateii, particularly if the 
donor potential is not Coulombic as is the case in the immediate vicinity of the donor nucleus. The consequence of 
this is that the donor electron binding energy given by this wave function is E = 28.95meV, significantly lower than 
the experimental value of E — 45.5meV 12,13 . This discrepancy is thought to arise from the deviation of the donor 
potential from a purely Coulombic potential as well as the effect of a non-static dielectric constant in silicon 14 . It is 
expected however that at distances of more than approximately an effective Bohr radius from the donor nucleus, the 
donor potential should be approximately Coulombic and so Eq. [5] will provide a good description of the true donor 
electron in this region 8 . Thus the Kohn-Luttinger form of the donor-electron wavefunction should be adequate for 
the purposes of calculating exchange energies for donor separations in the range considered in this article. 

A plot of the donor electron wave function along three directions of high symmetry, calculated using effective Bohr 
radii as determined by Koiller et aim, is shown in Fig.|2for a donor placed at a substitutional lattice site. The coefficients 
^4k M ,G, were calculated using the simple local empirical pseudopotential method as outlined in reference 1 ^, and a basis 
of 125 states. This method accurately reproduces the electronic band structure for silicon, particularly in the region 
of interest for this calculation, the lowest conduction bands, as obtained using more complicated non-local pseudo- 
potential techniques. The donor electron wave functions obtained in this manner clearly display oscillations produced 
by the interference between the Bloch functions located at the different conduction band minima. The wavefunctions 
are real, and in the [111] direction slightly asymmetric, the asymmetry being a consequence of the presence of the 
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FIG. 2: The solid line shows Kohn-Luttinger wave function for a phosphorus donor electron in silicon, plotted along directions 
of high symmetry within the crystal. The dotted line shows an isotropic Is hydrogenic wave function, with an Bohr radius of 
20.13 A. 

second sub-lattice. In the [111] orientation the silicon atoms are nor evenly spaced, and so the neighbouring silicon 
atom on one side of the phosphorus donor is much closer than that on the other side. Superimposed over the donor 
wave function is an isotropic Is hydrogenic envelope with a Bohr radius of 20.13 A , reflecting the effect of the 
superposition of the non- isotropic envelope functions F^. 

III. THE HEITLER LONDON FORMALISM 

The two electron Hamiltonian for a system of two donors separated by a vector R, in effective Rydberg units is: 



H 



2m 2 

„2 
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2m 2 
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4?re n - r 2 



(6) 



In the standard Heitler-London approximation^ one assumes that the lowest energy two electron wavefunction of the 
two-donor system is simply a correctly symmetrised superposition of single electron wavefunctions centred around 
each of the two donors 



1 



tf±(ri,r 2 ) = -= (#(ri - R/2)*(r 2 + R/2) ± *(n + R/2)*(r 2 - R/2)) , 
v2 



(7) 



where the two donors are located at positions ±R/2. This approximation is asymptotically exact and should hold 
for separations greater than the effective Bohr radii, |R| >> a±,au. The anti-symmetry of the fermion wavefunction 
then tells us that the exchange splitting, the difference in energies between the spin singlet and triplet spin states, is 
simply equal to the difference in energy between the states that is -Etripiet — ^smgict = E- — E + . 
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FIG. 3: Calculated exchange couplings as a function of donor separations along three high symmetry directions. In each case 
we plot with a solid line the results obtained when the periodic part of the Bloch function, u(r), is included, a dashed line 
indicates the results obtained with u(r) = 1. In each plot these lines cannot be separately resolved. The dotted line gives the 
results obtained when using the approximations of KHD and the dotted-dashed line gives the exchange coupling calculated 
using a Is hydrogenic orbital with an effective Bohr radius of 20.13A. The * denote the positions of lattice sites. 



In this article we present our results in terms of the exchange, or J coupling, in the exchange term of the effective 
spin Hamiltonian, Jcrf .erf, which is common to the quantum computing literature. We make this decision despite 
the observation that it is the exchange splitting, the energy difference between the single and triplet states, that is 
most commonly presented in the solid-state literature. These quantities are related by the expression J = (-©triplet — 

^singlet )/4. 

In the Heitler-London formalism the exchange coupling can be expressed in terms of matrix elements of the Hamil- 
tonian one can rewrite this as: 

J(R) = ±S(R) 3 (Ii(R) - J 2 (R))/(1 - S(R) 4 ), (8) 
where the overlap integrals are given explicitly by: 

Ji(R) = y > d 3 r 1 d 3 r 2 ^*(r 1 +R/2)**(r 2 -R/2)iJ*(r 1 +R/2)1'(r2-R/2), 

J 2 (R) = yd 3 rid 3 r 2 **(ri+R/2)**(r 2 -R/2)i/*(r 2 +R/2)*(ri-R/2), 

5(R) = y*d 3 r*(r + R/2)**(r-R/2). (9) 

Computation of these expressions using the Kohn-Luttinger wave function without approximation is rather tedious, 
given the large number reciprocal lattice vectors it is necessary to sum over to evaluate the periodic part of the Bloch 
functions for each of the six Bloch states in the donor electron wavefunction. We have used an adaptive Monte-Carlo 
quadrature program to evaluate all overlap integrals, taking due care to obtain reasonable precision. 

In Fig. |3 we plot the exchange energy as a function of donor separation in each of the major lattice directions. 
For comparison the result using an isotropic hydrogenic wave function is also given. We compare our results with 
those obtained using the method of KHD in their original papesi, who used two major approximations in order to 
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make the calculation more tractable, some of which are discussed in detail in reference a 2 i 17 . First, the Heitler-London 
expression for the exchange energy is approximated by the Coulomb exchange integral alone, an approximation that is 
asymptotically quite good, however it is this approximation that is responsible for the difference between our result and 
that obtained by KHD in referencei. In reference^, the authors have used the complete form for the Heitler-London 
expression, and obtained results that match those presented here. 

The other approximation made is to effectively ignore the contribution of the periodic part of the Bloch functions 
it M (r) = 1; this is an excellent approximation as can be seen from the figures in which it is almost impossible to 
distinguish between the results obtained in this approximation and those for which the detailed Bloch structure was 
included. The ability to ignore the periodic part of the Bloch structure in the computation of the Heitler-London 
integrals significantly reduces the complexity of the calculation by eliminating the sums over reciprocal lattice vectors. 
It is worth examining this approximation in more detail as we find it is effectively the same approximation as is made 
in effective mass formalism to obtain the Kohn-Luttinger form for the donor electron wavefunction. 

The Heitler-London integrals are all of the form 

I = J **(ri - R/2)tf*(r a + R/2)V(n, r 2 )*(r 2 - R/2)tf (n + H/2)dv 1 dr 2 
= f J2 X! Fi(ri-B./2)F j {T2+R/2)Fi(T2-R/2)F m (T l + R/2)V(Ti,T 2 ) 

x e i(k«+k m -kj-k,).R/2 e i(k m -k«).r le i(k,-k i ).r 2 

x A^^A^^^^A^G^^^^e^^^^dv^. (10) 

The various potentials represented in the above expression as V(ri, r 2 ) are all impurity potentials of a similar nature 
to those encountered in the effective mass formalism, and so we expect that they are sufficiently slowly varying that we 
can safely ignore the rapidly oscillating terms in the above integrand, in exactly the same manner. This immediately 
gives 

1 = J E F '( r i - R/2)^(r 2 +R/2)^(r 2 - R/2)i=)(r x + R/2)^(r x , v 2 )S M ^-^l 2 dr x dr 2 , (11) 

id 

which is equivalent to setting (r) = 1 in Eq. The "rapidly oscillating" terms include the exponentials containing 
the values of k at the conduction band minima, which are separated by a minimum magnitude of <5k = 2\/2 x 
0.857r/d. This separation is sufficiently large to allow us to ignore terms except those for which i — m,l = j. Since 
Sg^Ic G^k g = 1' this auows us to ignore the periodic part of the Bloch functions in the Heitler-London Integrals. 



IV. FABRICATION 



The observed extreme sensitivity of the exchange coupling on the relative orientation of the two phosphorus donors 
sets stringent requirements on the placement of donors for any quantum computer architecture that is reliant on the 
exchange interaction to couple qubits. In this section we discuss these implications for the construction of a Kane 
nuclear spin quantum computer. 

Fabrication of the Kane solid-state quantum computer is being pursued along two parallel pathsi^. In the so called 
"bottom up" approach, individual phosphorus donors are effectively placed with atomic precision on a phosphorus 
surface by application of phosphine gas to a hydrogen terminated silicon surface in which individual hydrogen atoms 
have been removed at the proposed donor site using a scanning-tunnelling microscope tip. The hydrogen mono-layer 
is then removed and the surface overgrown with phosphorus 19 . Although this approach has not been fully developed, 
it is likely that it will be able to produce an array of donors located at precise fee substitutional sites, to within a few 
lattice spaces. Small displacements are inevitable however and can have dramatic effects on the exchange coupling 
between the donors. This is illustrated in Fig. 0] where we show a plot of the exchange coupling as a function of the 
magnitude of the displacement \8\ of the second phosphorus dopant from its desired position at a fee substitutional 
site a distance of 200. 9lA, in the [100] direction, from the first dopant. These dopants have been displaced by no 
more than four fundamental lattice spacings on either of the two fee lattices that make up the silicon matrix, however 
the exchange coupling varies by more than an order of magnitude emphasising the need for precise placement. 

The second approach, known as the the top down approach calls for the implantation if the phosphorus donors into 
the silicon using an ion beam of phosphorus ions incident on the substrate^. The precision of such a technique in 
the placement of dopants is fundamentally limited by the phenomenon known as "straggling" , whereby the implanted 
ions scatter from the nuclei of the host silicon atoms. Simulations of this process for a beam of 15keV phosphorus 
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FIG. 4: Calculated exchange couplings for donors at fee lattice sites that are displaced by a vector 8 from their ideal separation of 
200. 9lA in the [100] direction. The couplings are plotted as a fraction of the expected exchange coupling J(200.9lA) = 0.18^ieV. 
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FIG. 5: A plot of the logarithm of the integrated probability for the exchange coupling between donor electrons for donors im- 
planted 200A apart in the [100] direction, using a 15keV phosphorus ion beam. The donors are assumed to take a substitutional 
fee lattice site. 

ions implanted into silicon, an energy appropriate for an implantation depth of approximately 200A into the silicon 
substrate, give a roughly Gaussian distribution for the final position of the dopant with a variance of about 90A in 
the transverse direction, and 100A in the longitudinal direction. Using these data, for two dopants implanted 200A 
apart in the [100] direction, we have calculated the distribution of exchange couplings between the donor electrons, 
assuming that after thermal annealing the phosphorus donors take up a position on a fee substitutional site. Fig. [5] 
shows a plot the integrated probability, JJ° P(J')dJ' , that is the probability that the exchange coupling greater than 
a value J, for two dopants implanted in the manner just described. 

In the case of electron spin quantum computers where it is the exchange energy that directly couples neighbouring 
qubits, a coupling of 50/ieV corresponds to a gate time of the order of 10GHz. The top down approach can in 
this case produce qubits coupled by extremely fast gates with a very high probability. In nuclear spin quantum 
computers however the exchange interaction only mediates the inter-qubit coupling and the resulting gates are slower. 
The original Kane proposal calls for an adiabatic implementation of a controlled-NOT gate that requires a exchange 
coupling of approximately 50/ieV^Sl. Additional proposals exist for non-adiabatic implementations 2 ^^ for which 
controlled-NOT gates can be accurately implemented with a lower exchange coupling. Fig. [S] shows the probability 
that the bare exchange coupling for implanted donors being greater than 50/ieV, is about 6%. What is more important 
than the bare coupling however is the values of the exchange coupling that can be achieved with the application of 
J-gate biases. 

We note that KHD have performed calculations that suggest a method for overcoming the strong dependence of the 
the exchange coupling on the donor orientation 2 . They calculate the exchange coupling, in the same approximations 
as discussed earlier, for phosphorus donor electrons in uni-axially strained silicon. The strain is a product of the the 
silicon host being fabricated on a layer of Sii-^Gea, and is found to suppress the oscillations in the coupling for donors 
that lie within a plane perpendicular to the direction of the uni-axial strain. The coupling remains highly sensitive 
to displacements away from this plane. 
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FIG. 6: Potential produced inside the device by a voltage bias of IV applied to a metallic J-gate. Here x is the inter-donor 
axis, the donors are located at x — ±100A and at z — 0. The potential is independent of the ^-coordinate, mimicking a very 
long electrode. The dashed line is the one-dimensional potential used by Fang et al with /i = 0.15 and shifted to match our 
potential at x, z — 0. 

V. EXTERNAL CONTROL OF THE EXCHANGE COUPLING 

Inherent in the Kane proposal for a solid state quantum computer is the ability to control the exchange coupling 
of neighbouring donor electrons through the application of voltage biases to a control J-gate placed on the substrate 
surface between the position of the phosphorus donors, as illustrated in Fig^ In this section we calculate the effect on 
the exchange coupling on the J-gate bias. Some results have been reported on similar calculations. Fang et atL used an 
unrestricted Hartree-Fock method, which avoids some of the approximations inherent in the Heitler-London approach, 
to calculate the exchange coupling between phosphorus donor electrons as a function of a simplified electric potential. 
They use a trial wavefunction of the form \&(r) = F(r)0k fl (r), including only one of the six degenerate conduction 
band minima, with a spherical envelope function, that is one of the same form as Eq. but with a± = an = 20A. As 
a result of their single minimum approximation they cannot possibly reproduce the oscillatory nature of the coupling 
which results from interference between the terms localised at the six minima. They model the potential produced 
inside the device by the control gate as one- dimensional potential of the form Vj{x) — (i ^"^j^j^^ 2 ^ x 30meV, 
where x is the distance is along the direction between the two donors which are situated at x = ±i?/2, as defined in 

Fig.m 

In our calculation we compute the potential produced inside the device due to the application of a voltage bias to a 
metallic J-gate by numerical solution of the Poisson equation using a commercial package, TCAD 24 , designed for the 
semi-conductor industry. In Fig. we plot the potential parallel to the inter-donor axis, for several values of the z 
co-ordinate, which denotes the distance below the oxide barrier, located at z = — 200A. For more information on the 
details of the potential calculation see ReS^. We compare the potential obtained for a J-gate voltage of IV, with the 
one-dimensional potential of Fang et al and see that the two potentials agree well for z = 0, the plane of the donors, 
provided we set fi = 0.15. However we see that for this value of \x the two potentials are quite different in the planes 
50A above and below the donors. We must therefore conclude that the one-dimensional approximation is not a good 
one. 

To calculate the effect of the applied potential on the electron we expand the wavefunction as follows: 

*(r-R/2,Vj) = ]T c„, i , m (K/)^F;^ m (r-R/2)e 4k -( r - R / 2 ), (12) 

n,l,m M 

where the F™' l ' m (r — R/2) are non-isotropic hydrogenic envelope functions defined by Faulkner—. The n is a label for 
the six degenerate conduction band minima and determines the direction of the anisotropy of the envelope function, 
for example F± z ' m (r) = F n,l,m (x,y,'jz), is the hydrogenic function with a Bohr radius a — a±, and 7 = an/a_\_- 
The orthonormality of this basis is enforced by the e l ( k A»- k ^) ( I '- R / 2 ) terms which appear in the matrix elements and 
due to their rapidly oscillating nature average to zero unless k M = k„. The coefficients c ni i_ m arc determined by 
direct numerical diagonalisation of the Hamiltonian in the presence of the electrostatic potential. We use a basis 
of 140 states, which includes all states with principle quantum number up to and including n = 7, and we re-scale 
the energies such that the unperturbed ground-state energy gives the experimentally observed value of 45.5meV for 
phosphorus donors in silicon. 
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FIG. 7: Exchange coupling as a function of donor separation for donors oriented along the [100] direction and various J-gate 
biases. 
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FIG. 8: Exchange coupling as a function of J-gate bias for donors at fee lattice sites separated by 200. 9lA in the [100] direction. 



The donor electron wavefunctions obtained in this way are then used in the Heitler-London formalism to calculate 
the electron exchange energy, as a function of both the donor separation, and applied J-gate bias. The results are 
plotted in Fig. for donor separations from 80-120A in the [100] direction. We see that application of a positive 
bias will tend to increase the exchange coupling, while a negative bias decreases the strentgh of the coupling. It is 
also worth noting that the relative increase in the exchange coupling produced by a given bias increases with the 
donor separation. Also the change in coupling is strongest at the peaks, the application of a positive bias enhances 
the oscillations in the exchange coupling whereas a negative bias of magnitude IV inverts the oscillations such that 
points that were originally at peaks become troughs. 

The range of separations shown in these plots is far less than the approximately 200A donor separation called for 
in the original Kane proposal- we have plotted values over range for ease of comparison with our previous results, 
and more particularly those of KHD. In Fig. IHlwe plot the exchange coupling as a function of the J-gate bias for 
donors separated by 200. 9lA in the [100] direction. At this separation a IV J-gate bias can increase the the exchange 
coupling by over two orders of magnitude up to a value of approximately 30/ieV. The exchange coupling can also be 
reduced via application of a negative bias. 

Although it is difficult to compare our results with those of Fang et al due to the different potentials, the shift of 
an order of magnitude obtained for separation of 200 A and a bias of IV is greater than the shift they predict. For 
a value of /i = 0.15 Fang et al predict a shift of less than an order of magnitude. This discrepancy can in part be 
attributed to the 1 — D approximation for the gate potential which does not account for the fact that the magnitude 
of the potential increases as the J-gate is approached from below (in the z-direction). We note that it may not be 
possible to apply voltages much greater than one volt to these nano-scale devices without exceeding the breakdown 
fields of the various materials involved 25 . 
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VI. DISCUSSION AND CONCLUSIONS 



In agreement with the approximations implicit in the calculations of donor electron exchange energy of Koiller et 
aim we show explicitly that the periodic part of the Bloch functions in the donor electron wavefunctions can be ignored 
in the Heitler-London integrals, greatly reducing the complexity of the calculation. This has allowed the calculation of 
the intrinsic exchange coupling for a large number of donor pairs distributed according to the probability distribution 
for donors implanted 200A apart by 15keV ion beams. In this manner we have determined the probability distribution 
of the exchange coupling for the top-down approach for the fabrication of a Kane solid-state quantum computer. 

In addition we have investigated the application of control gate biases to increase the exchange coupling between 
donor pairs using both realistic potentials and realistic donor electron wavefunctions . We find that as expected the 
exchange coupling can be increased by the application of a positive bias, and decreased with a negative bias. Over 
the range of donor separations investigated it was found that the relative increase of the exchange coupling produced 
by a given bias increases with the separation, and for donors separated by 200A in the [100] direction a IV bias can 
increase the coupling by over two orders of magnitude. 
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